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One of the main goals of relativistic heavy ion collisions is the determina- 
tion of the nuclear equation of state. At high energies, semiclassical cascade 
models in terms of scattering hadrons have proven to be rather accurate in ex- 
plaining experimental data. Therefore it is of fundamental interest to extract 

the equation of state from such a 
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Figure 1: 'Equation of state' of infinite nu- 
clear matter, calculated with URQMD (pre- 
liminary). Shown is the energy density as 
a function of temperature (extracted with 
a least square fit to the baryon momen- 
tum spectrum) at fixed net-baryon density of 
p = 0.16 fm — . The equations of state of a 
Hagedorn-gas, of a quark-gluon plasma, and 
of an ideal gas of nucleons and ultrarelativis- 
tic pions are also depicted. 



microscopic model, i. e. to investi- 
gate the equilibrium limits and bulk 
properties, which are not an explicit 
input to the non-equilibrium trans- 
port approach with its complicated 
collision term (unlike e. g. in hy- 
drodynamics). In Fig. [j] we study 
the thermodynamic properties of in- 
finite nuclear matter with the Ultra- 
relativistic Quantum Molecular Dy- 
namics 1.0/3 (URQMD), a semiclas- 
sical transport modeEI. The model 
is based on classical propagation of 
hadrons and stochastic scattering (s 
channel excitation of baryonic and 
mesonic resonances/strings, t chan- 
nel excitation, deexcitation and de- 
cay). For this study the potentials 
have been switched off. To sim- 
ulate infinite hadronic systems we 
construct a box of volume 250 fm 3 
with periodic boundary conditions, 
and initialize 40 nucleons — that is 
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Figure 2: Transverse momentum spectrum of directly produced photons in Pb+Pb collisions 
at 160 GeV/u calculated with URQMD (left). The contributions of the different processes 
are shown. The resulting spectrum is compared with hydrodynamical calculations (right). 
For the photon source, an equation of state with and without a phase transition is assumed. 



ground-state nuclear density — randomly in phase space, varying the total 
energy density. After the system has equilibrated according to the URQMD 
simulation, the temperature is extracted by fitting the particles' momentum 
spectra. Alternatively, the temperature can be extracted from the relative 
abundances of different hadrons, e. g. the A/N ratio, which should yield the 
same temperature. This has been checked in the current simulation. The result 
of this procedure is plotted in Fig. |[ It appears that the energy density rises 
faster than T 4 at high temperatures of T rj 200 — 300 MeV. This indicates 
an increase in the number of degrees of freedom. It may be interpreted as 
a consequence of the numerous high mass resonances and string excitations, 
which in a way release constituent quark degrees of freedom (but, of course, no 
free current quarks as in an ideal QGP). Investigations of equilibration times 
and relative particle and cluster abundances are in progress. Moreover, the 
admittedly poor statistics have to be improved. 

Experimentally, one can access the equation of state of strongly interacting 
matter by measuring thermal electromagnetic radiation emitted in heavy ion 
collisions. We have calculated direct photon production in Pb+Pb collisions 
at 160 GeV/u within the framework of URQMD 1.0 (3. To improve statis- 
tics, the total meson-meson cross sections were fixed to 15 mb (independent of 
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fJs) . For the sake of comparison with earlier hydrodynamical calculations (see 
□ and refs. therein) we considered only the .processes irp — ► 777, irrj — ► ttj and 
7T7T — > ^7. The amplitudes are taken fromtl As can be seen in Fig. || (left) 
the 7r + p — > 7r + 7 process dominates the spectrum in the transverse momen- 
tum range considered. Fig. || (right) shows that the direct photon slope from 
the microscopic calculation equals that from the hydrodynamical calculation 
without a phase transition in the equation of state of the photon source. In 
the case of a phase transition the slope becomes significantly steeper due to a 
lower temperature at given energy density. The lower absolute yields of the 
direct photons from the URQMD may indicate that the meson abundance at 
early times is smaller than in the hydrodynamical model of 13, although the 
final yields are similar (due to higher mass mesons and baryon resonances). 
Thermal and in particular chemical equilibrium may not be established within 
URQMD in the early stage of the collision (where the large transverse momen- 
tum photons are produced). In contrast, in the hydrodynamical calculation 
local thermal and chemical equilibrium of the produced particles is assumed 
from the |£ery beginning. Experimental studies on direct photons, under way 
at CERNo, will help to settle these questions in the near future. 
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